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Abstract 

A stratified charge engine employing the Texaco Controlled 
Combustion System has been operated over a large range of load 
conditions on iso-octane, methanol, and a wide boiling point 
(100-600°F) residual fuel. Basic performance, emissions and 
combustion parameters were measured over a range of overall 
equivalence ratios from ^=0,1 -1.0 at three engine speeds; 
1500, 2000 and 2500 RPM. The basic performance and emissions 
data were found to vary little between iso-octane and residual 
fuels, and to compare very well with similar data collected on 
the same engine design at other research facilities. The engine 
operation on methanol was not entirely satisfactory due to an 
improper match between the specific fuel injection system used 
for these experiments and the design requirements imposed by the 
much lower heating value and higher stoichiometric fuel-air 
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mass ratio of methanol. 

A direct, online data acquisition system, based on a 
Digital PDP 11/10 computer was developed to obtain accurate 
pressure-crankangle data for further combustion and thermo- 
dynamic studies. The acquisition program also computes a mean 
pressure - crankangle diagram and the statistics associated 
with cycle to cycle variation. The mean pressure - crankangle 
data is then integrated to compute indicated mean effective 
pressure. The opportunity to analyze pressure - crankangle data 
in this way substantially improves the accuracy and speed of 
data collection. 

A simple thermodynamic model based on homogeneous charge 
engine combustion has been modified to compute the heat release 
and fuel fraction burnt from the pressure - crankangle data. 

The problems associated with calculation of these parameters in 
diesel or stratified charge engines are discussed. Recommendations 
are made for further development of the online data acquisition 
system and the thermodynamic model. 
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I INTRODUCTION 



A stratified charge engine is defined as a spark ignition, 

internal combustion engine with a non-uniform fuel air mixture 

in the combustion chamber. Stratified Charge engines have been 

recognized for good fuel economy potential, low emissions and 

( 1 )* 

the ability to bum a wide range of fuel types. In this 

thesis, we analyse the performance of an engine based on the 
Texaco Controlled Combustion System (TCCS) , 

Several features specific to this design make the engine 
a strong candidate for small, low power applications. Load 
control can be achieved by changing the amount of fuel in- 
jected and inlet air throttling is unnecessary in most ap- 
plications. As a result, the engine can be run at very lean 
overall equivalence ratios giving excellent fuel economy and 
low emissions. Fuel is injected late in the compression stroke 
just before the combustion process and is ignited with a spark 
discharge. The residence time at elevated temperatures is 
therefore shorter than the time for compression ignition and 
hence the engine does not display either octane or cetane require- 
ments. 

As illustrated schematically in Figure 1, the Texaco 
Controlled Combustion System uses an open combustion chamber with 
high air swirl, direct fuel injection and electronic ignition. 

Air swirl is generated by the inlet air flow, and amplified 

^Numbers in parenthesis refer to the bibliography at the end of 
this paper. 
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( 2 ) 



during the compression stroke by constraining the vortex in 
the combustion chamber. The combustion chamber is essentially 
a cup with a cylindrical upper section and a toroidal bottom 
formed in the head of the piston. Fuel is injected with a 
Roosa Master Pencil Nozzle with a flat seat and a single hole 
orifice as shown in Figure 2, The positive ignition system 
uses a high energy multiple spark unit with controlled duration. 

The spark plug electrodes are carefully aligned to promote the 
formation of a stable flame front. 

High pressure injection of fuel into the swirling air 
begins near the end of the compression stroke. Air swirl in 
the combustion chamber promotes mixing and controls the pene- 
tration and trajectory of the fuel spray in the cup. The com- 
bustible fuel-air mixture formed by turbulent mixing and air 
entrainment in the fuel jet is then ignited by the long duration 
spark discharge and burns downstream of the spark plug as shown 
in Figure 3. 

(3) 

Following B.C, Jain we divide the combustion process 
into three stages; a rapid combustion phase controlled by the in- 
jection rate and a slower "burn up" phase which is controlled by 
the rate on air entrainment and mixing of burning products and 
air downstream of the spark plug# and a heat transfer dominated . phase 
which follows after mixing is complete.' This sequence is illustrated 

on the diagram shown in Figure 4. During an isentropi£ compres- 

v 

sion or expansion of an ideal gas the value of PV remains 









p 




































- 12 - 



constant. After a delay covering the jet transit time for the 
injector to the flame front, the value of PV rises rapidly. 

This rapid combustion phase appears to be controlled by the in- 
jection rate. After the last fuel injected passes the spark 
plug, the rate of change of PV is substantially slower and the 
rise is controlled by the rate of mixing of the plume of rich 
products with the surrounding air and residual gas. This 
phase ends when all mixing is complete or the exhaust valve 
opens. Any fall in the PV curve prior to the exhaust valve 
opening can be attributed to heat losses. 

This research project is a part of a larger program which 
includes work on a jet mixing model, a performance model and 

( 4 ) 

photographic studies with a rapid compression machine. 

The following areas of research are covered in this thesis; 

(1) The completion of the engine test setup and the 
development of all necessary instrumentation to 
record the variables of interest. 

(2) TCCS engine performance and emissions are pre- 
sented for the three test fuels over a wide 
range of operating conditions. Differences in 
fuel characteristics are also presented. 

(3) The problems associated with heat release cal- 
culations in stratified charge and diesel engines 
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reviewed, with the further development of an 
existing thermodynamic model to predict heat 
release and mass fraction of fuel burnt outlined, 

(4) Log P vs Log V diagrams and plots summarizing the 
output of the thermodynamic model are presented 
for a matrix of comparable test data, 

(5) The development of computer programs to process 
raw data and accomplish online pressure data 
acquisition, along with listings of all computer 
programs are included as appendixes. 



-14- 



II TEST ENGINE AND INSTRUMENTATION 



Single Cylinder Test Engine 

The engine used in these experiments is located in the 
Sloan Automotive Laboratory at MIT, and is arranged as shown 
in Figure 5 . The single cylinder test engine is based on a 

CFR - 48 crankcase that has been modified to accept a cylinder 
sleeve assembly, head, piston, crankshaft and overhead cam and 
valve train assembly for the 3 7/8 inch bore by 3 7/8 inch 
stroke LIS - 183 TCCS geometry, as shown in Figure 6 . 

Engine specifications and dimensions are shown in Table 1. The 
engine is coupled to a dynamatic eddy current dynamometer 
equipped with a hydraulic scale, as shown in Figure 7 . The 

basic instrumentation is listed in Table 3 and whenever possible 
redundant measurements have been introduced to provide alternate 
data sources and to qualify experimental results. 

The basic engine support facilities were constructed by 
Lazarewicz ^ ^and follow standard practices as shown in Figures 8 
through 11 . The engine cooling system is arranged to enable 
heat rejection measurements. A rotameter and throttling valve 
on the return line from the engine allows both flow regulation 
and measurement while holding system pressure above 5 psig, as 
shown in Figure 8 . Maximum water temperature at the engine cy- 
linder outlet was held at 190 ±5°F. The lubrication 
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system shown in Figure 9 consists of two separate loops; a low 
pressure circulating loop for temperature control and a high 
pressure bearing feed and filter loop. The operating oil temp- 
erature and pressure were held at 165°F and 42 psig respectively. 

A pressure alarm incorporated in the filter system activates a 
siren and cuts off the fuel supply and ignition system when 
pressure falls below 20 psig. Inlet air flow is measured with 
an ASME square edge orifice with flange taps and water manometers 
as shown in Figure 10 . Following recommendations made by 
Lazarewicz, the inlet and exhaust systems were rebuilt and the 
inlet settling tank and air heater were closely coupled to the 
engine with a short inlet pipe. 

Injection and Ignition Systems : 

Fuel system specifications appear in Table 3 and fuel flow 
is measured gravimetrically as shown in Figure 11 . Transfer 
pump pressure is held at 25 psig # Fuel injector leakoff and 
the injector pump bleed are returned to the fuel reservoir mounted 
on a standard laboratory scale. 

Ignition and injection timing, intensity and duration 

are monitored with a 565 Techtronix Oscilloscope. Ignition timing 
can be precisely set by a means of a vernier scale and adjustment 
arm. Injection timing is varied through the use of an American 
Bosch TMB - 12 Manual Timing mechanism. A pressure transducer 



-16- 



is mounted just ahead of the fuel injector in the high pressure 
supply line and the Needle Lift Indicator shown in Figure 12 
has been developed to replace the standard fuel injector needle 
lift and cracking pressure adjustment assembly^ Both 

outputs are displayed on the oscilloscope cathode ray tube 
along with the cylinder pressure crankangle markers and BDC 
reference pulse as indicated in Figure 13 . As presently de- 
signed, resolution of fuel injector timing is limited to 2 CA° 
by the time base necessary to include one entire revolution on 
the CRT display. The first peak in the pressure trace cor- 
responds to the point of initial needle lift, however; as in- 
dicated by the needle lift trace, the injector does not open ap- 
preciably for about 2 CA°. Determination of injector cutoff and 
any secondary injection is only possible with the Needle Lift 

i 

Indicator. Substantial signal processing problems occur with the 
linear displacement amplifier and a new design is 
recommended. 

Emission Sampling System : 

Engine emissions are measured using the Sloan Laboratory 
Exhaust Gas Analysis Cart (Table 3). The exhaust sample is re- 
moved from an engine exhaust tank consisting of several cylinder 
volumes and pulled through a heated teflon line and filter to the 
Gas Analysis Cart. Stainless steel pipe is used between the 
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engine and exhaust surge tank to reduce the potential for 
hydrocarbon reactions induced by M rust ,f and the elevated 
temperatures of the exhaust flow. An additional spun glass 
particle filter was installed in the sampling line when it was 
discovered that at high loads and equivalence ratios near <f> = 1.0, 
the gas analysis equipment was severely contaminated with carbon. 
This technique appears to have solved this immediate problem. 

The development of a computer program to calculate basic 
engine performance and emissions from the recorded data is 
described in Appendix I with the complete program listing in 
Appendix III. 

Pressure Volume Measurement : 

Accurate combustion pressure and volume measurements are 
absolutely required for mathematical engine simulation, cal- 
culation of heat release rates, engine pumping losses and the 
comp ilia ticn of statistics associated with cyclic variations and 
peak pressures. Acceptable pressure volume records can only be 
obtained when the entire monitoring system receives careful 
attention. First, high resolution signal recording equipment is re- 
quired if the effort expended to obtain accurate pressure and volume 
measurements is to be worthwhile. Oscilloscopes were used in these 
experiments for qualitative analysis but they lack the accuracy 
required for quantitative resolution. The real time digital 
computer with analog to digital converters can provide the 
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resolution necessary; and for these experiments, an online 
digital data acquisition was developed for Sloan Laboratory 
PDP 11 Computer Facility, The analog to digital converter 
provides resolution to within 0.48 psi over a range of 0-1000 psia, 
well in excess of current engine pressure transducer performance 
as described in following paragraphs. The sampling intervals 
were 5 CA°. Appendix II provides details of the system hard- 
ware and the development of the software used for online 
computer sampling. Appendix III contains a listing of the 
pressure-crankangle data management program and the assembly 
language program that actually performs the data acquisition. 

In these experiments the cylinder head geometry prevented 
the use of a large diaphragm, water cooled pressure transducer 
and a Kistler 609A piezoelectric pressure transducer was chosen. 
Considerable experimental art is required to obtain acceptable 
performance from available pressure measuring equipment, in- 
cluding this specific unit. For example, quartz piezoelectric 
transducers are high impedance devices and contamination of the 
electrical connectors can significantly degrade performance. All 
connections must be thoroughly cleaned with a freon base solvent 
and sealed with heat shrink tubing. Transducers are also subject 
to the thermal cycling that is fundamental to the combustion 
process in engines. In a chopped flame test by Jain and 
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Lazarewicz^\ the 609A transducer showed an apparent 6 psi 
response when directly exposed to an acetylene flame at typical 
engine frequencies. This response was reduced with a coated 
diaphragm and a thin coating of silicone rubber was applied as 
recommended in the literature^ ^ \ The transducer was in- 

stalled in a recessed adapter inserted through the water jacket 
of the cylinder head. The. adapter cavity is designed to 
minimize attenuation and protect the transducer from engine 
temperature fluctuations. The engine coolant also serves to cool 
the transducer. 

Careful preparation of the transducer is wasted if the cylinder 
volume is not known with similar accuracy. The cylinder volume 
is computed from engine dimensions and recorded crankangle data. 
Cylinder clearance volumes were determined by careful measure- 
ment of pertinent engine dimensions. The piston top dead center 
was located and the flywheel position pointer adjusted with a 
depth micrometer as recommended by Lancaster^ ^ . A crankshaft 
driven rotary pulse generator, supplying 720 pulses plus a marker 
every revolution was then aligned with the flywheel with an 
accuracy of approximately 1/4 degree. The alignment was tested 
by analyzing pressure crankangle diagrams and Log P - Log V plots 
of motoring runs as shown in Figures 14, 15 and 16 . Lancaster 

provides a detailed explanation and interpretation of the plot 
(7) 



orientations 
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A clamped disk balanced pressure indicator was also 
mounted through the engine water jacket to provide a technique 
for dynamic calibration for the cylinder pressure transducer. 

The balance pressure indicator is a pressure activated switch 
with a reference pressure applied to one side of a thin membrane 
disk and the cylinder pressure to the other. The balanced 
pressure indicator is connected to the oscilloscope display as 
shown in Figure 17 . When the cylinder pressure rises above the 
reference pressure plus the disk contact pressure, the disk is 
deflected to ground the center electrode. The change in potential 
is converted by the cathode ray tube grid modulator to momentary 
changes in signal intensity on the oscilloscope display as shown 
in Figure 13 . These pulses are then used to dynamically 
calibrate the signal from the piezoelectric transducer. 
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III BASIC PERFORMANCE AND EMISSIONS 

The multifuel capability of the TCCS engine was in- 
vestigated using methanol, iso-octane and a wide boiling 
point (100-600°F) fuel. Properties of these test fuels are 
summarized in Table 4. Performance and emissions data in 
Figures 18 through 39 represents engine operation over the range 
of fuel-air ratios and engine loads summarized in Table 5. The 
engine was naturally aspirated and exhausted to ambinet pres- 
sure throughout the test series. All performance data was 
measured with injection timing set for maximum brake torque 
and ignition system timing set to commence 2 CA° prior to the 
start of injection. The injection duration exceeded the ig- 
nition duration of 20 CA° except at light load. The overall 
equivalence ratio was used as the absissa in Figures 18 through 
38. The equivalence ratio was determined using two techniques; 

a) the measured fuel and air flow and b) calculated from the 

/ s \ 

exhaust gas composition using the method of Stivender . Data 
was considered reliable when the difference between these two 
values was less than 0.025. A deviation greater than this was 
always traced to operator error, air leaks or faulty equipment. 
The engine was not operated at fuel air ratios above stoichio- 
metric. Previous experience with this engine has shown that any 
performance gain above 4> = 1.0 are achieved with substantially 
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increased hydrocarbons, CO and degraded fuel economy 

Engine Performance with Iso-Octane and Wide Boiling 
Point Fuel ; 

The indicated mean effective pressure (IMEP) versus 
equivalence ratio for iso-octane and the 100-600, wide boiling 
point fuel are shown in Figures 18 and 19. The maximum IMEP 
is not developed, with either fuel, in the range of equivalence 
ratios shown; rather, the effective upper limits for engine 
operation are determined by a "smoke limit” near ^ = 1.0. In 
addition, there is a distinct flattening of the power curve near 
stoichiometric conditions. The secondary dependence on RPM 
exhibited by the IMEP curves can be traced to several factors. 

The volumetric efficiency increases with speed in the range 
tested as shown in Figure 39. PV plots also show a slight de- 
crease in heat transfer with increasing speed during the heat 
transfer dominated phase of combustion. As will be discusses in 
the following chapter, the burning angle appears to decrease 
slightly with increasing RPM, and this would also serve to in- 
crease the IMEP. The data for iso-octane and 100-600 split 

differently with speed. Similar results have been observed by 
(9) 

Texaco v ' . This difference cannot be satisfactorily explained 
with the basic data comparison presented in this thesis and ad- 
ditional study to resolve this potential conflict is recommended. 

The indicated specific fuel consumption (ISFC) is shown for 
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iso-octane in Figure 20 and for 100-600 in Figure 21. The 
data sets are of similar character with a clear minimum at an 
equivalence ratio near $ = 0.3 and a steep rise at leaner 
equivalence ratios. This rise in ISFC at lean fuel air ratios 
is accompanied by increased cyclic variations and incomplete 
combustion; and appears to be a characteristic of the fuel 
injection system used with this engine. 

The indicated thermal efficiency (*? ) provides the best 
comparison of the actual combustion process since it properly 
accounts for the different heating values of the three fuels 
used in these experiments. The indicated thermal efficiency is 
equal to the reciprocal of the ISFC x lower heating value and 
the indicated thermal efficiency is equivalent for engine 
operation on both fuels, with a miximum value of approximately 
50% reached near <f> = 0.3 as shown in Figures 22 and 23. 

The volumetric efficiency (*K) changes with load to re- 
flect changes in the quantity, composition and state of the 
residual gas in the combustion chamber as shown in Figures 22 
and 23. The effect of engine speed is shown in Figures 22, 23 
and 39. The effect of load is shown in brackets on Figure 39. 

The exhaust temperature data follows the same trends as 
described for the IMEP data as shown in Figures 24 and 25. 

The friction mean effective pressure (FMEP) versus RPM 
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for the single cylinder test engine is shown in Figure 39. This 
data is representative and the actual FMEP showed little variation 
throughout the test series. 

Emissions with Iso-Octane and Wide Boiling Point Fuel : 

The TCCS concept is designed to burn the fuel immediately 
after injection in a fuel rich, mixing controlled plume in 
order to achieve multifuel capability and low emissions. In the 
course of these experiments, it was consistently observed that 
emissions, particularly hydrocarbon and carbon monoxide , are 
more sensitive to small variations in engine operating conditions 
than the basic performance data. Careful system timing is re- 
quired to obtain satisfactory emissions levels and the convention 
adopted places the start of ignition immediately ahead of in- 
jection. However, if timing is adjusted so that fuel injection 
preceeds ignition higher IMEP levels are achieved with a cor- 
responding significant increase in hydrocarbon emissions at 
equivalence ratios of <£> 0.6. This system sensitivity is thought 
to be the source for data scatter at high load conditions and all 
lines were drawn using a least square regression analysis 
technique. 

The increased cyclic variations and degraded emissions 
observed at very low load test conditions are thought to be 
caused by an injection - ignition system phenomena. Two 
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potential mechanisms for cyclic variations are advanced. High 
speed movies of TCCS combustion in a Rapid Compression Machine 
show that the ignition arc discharge time is short compared 
to the arc cycle time^^. This ignition characteristic 
suggests that the initial fuel jet may pass the spark plug in 
the interval between ignition pulses and a stable flame front 
may not be formed in the leading edge of the fuel jet. This 
mechanism introduces the possibility of cycle variation in 
the initial stages of combustion. At light load conditions 
the injection duration is less than the ignition duration, and 
with small injection quantities, the unburnt fuel vapor that 
passes the electrodes before a stable flame front is formed can 
rapidly mix to a fuel air ratio below the limit of combustion. 
This mechanism may partially account for the degradation of 
hydrocarbon emissions at low load. 

A second possible mechanism is traced to fuel injector 
characteristics. As indicated in Figure 13, needle lift is not 
always crisp and there is often a 2 CA° period at the start of 
injection when the needle opens only a small amount. This 
starting transcient may have substantially influenced the 
initial stages of fuel jet formation. As shown by Jain, a low 
momentum jet would be swept outside of the electrode radius by 
air swirl. The percentage of fuel vapor missing the plug 
electrodes would increase for light load conditions with the 
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smaller fuel quantities required. As in the previous case, this 
phenomenon introduces a mechanism for cycle variation and 
hydrocarbon and CO formation in the combustion chamber. 

Carbon monoxide emissions are shown in Figure 26 for iso- 
octane and in Figure 27 for 100-600. The lowest CO levels are 
obtained at an equivalence ratio near <f> = 0.45. The CO levels 
observed with iso-octane are approximately 50% lower than observed 
with 100-600 at a given equivalence ratio. The mean minimum 
value observed with iso-octane 4 gr/ihp-hr or 3.5 gr/ihp-hr 
less than that observed with 100-600. 

The hydrocarbon emissions show a sharp increase at equiva- 
lence ratios less than 4> = 0.3 for both fuels, as seen in Figures 
28 and 29. This sharp rise is accompanied by increased cyclic 
variations which are attributed to the injection and ignition dynamic 
effects discussed previously. A small rise in HC emissions is also 
observed near stoichiometric equivalence ratios. 

Nitric oxide emissions are shown in Figures 30 and 31. The 
trends indicated by the data have the same characteristic shape 
as demonstrated in homogeneous charge spark ignition engines; 
however, the peak occurs near an overall equivalence ratio of 
<t> = 0.6 whereas in a homogeneous charge engine the measured levels 
are generally higher and the maximum level occurs near stoichio- 



metric fuel-air ratios. 
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Methanol Performance and Emissions : 

Methanol was chosen as the third test fuel because it 
represented a severe test of TCCS multifuel capability. The 
fuel properties of methanol are significantly different from 
iso-octane or the wide boiling point fuel. Methanol has a 
much lower specific heating value and a higher stoichiometric 
air-fuel ratio and thus requires injection quantities nearly 
twice as large to achieve the same equivalence ratio. A 
Bosch injection system was used in these tests and the pump and 
nozzle geometry were selected for the reference fuels. No 
attempt was made to modify the pump or nozzle for the methanol 
experiments. As a result, the fuel system was operated off 
design. 

The performance and emissions data for methanol are pre- 
sented in Figures 32 through 38. Since the fuel system was not 
optimized for this fuel, wide scatter in emissions data was 
observed. The basic performance plots in Figures 32 through 35 
show the trends to be expected with a properly matched fuel in- 
jection system. Figure 32 showing IMEP data exhibits no 
sensitivity to RPM; however, "misfire" increased with higher 
RPM at high load conditions. The indicated specific fuel con- 
sumption, volumetric efficiency, and thermal efficiency data, as 
shown in Figures 33 and 34 has trends very similar to the 
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corresponding trends observed with oso-octane and 100-600 fuel. 
The indicated thermal efficency ) is almost identical to 
Figures 22 and 23; the maximum value of approximately 50% occurs 
at an equivalence ratio near <j) = 0.3. This demonstrates that 
the Texaco controlled combustion process is compatible with 
alcohol fuels; however, fuel system modifications are necessary 
to properly match the engine with this fuel type. 

Comparison of M.I.T. and Texaco Data 

The single cylinder engine emissions and performance data 

compare favorably with the TCCS engine data observed at the 

(9) 

Texaco Engine Development Laboratory . The maximum IMEP 
observed at MIT with 100-600 fuel exceeded the values observed 
at Texaco by nearly 5% at a given equivalence ratio. The two 
engines exhibited identical values of volumetric efficiency at 
each given operating point. Comparative plots by Lazarewicz 
showed good agreement with all emissions data except hydrocarbon 

emissions^. In these experiments, the level of hydrocarbon 

I 

emissions observed was considerably reduced when compared to 
values obtained by Lazarewics; however, the level is still higher 
than observed at Texaco. This difference in HC levels can be 
explained by the different sampling techniques used at the two 
facilities. The Texaco data was acquired in bag samples before 
analysis, while at the Sloan Laboratory a heated teflon line is 
used to sample directly from the exhaust tank. The heated 
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sampling line used at MIT eliminates the potential for con- 
densation of hydrocarbons, and in general, higher hydro carbon 
levels are measured in experiments with a heated sampling line. 
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IV ANALYSIS OF PRESSURE VOLUME DATA 



Accurate records of cylinder pressure-volume data are an 
important tool for evaluating performance. The indicated mean 
effective pressure and pumping mean effective pressure can be 
computed by integrating the average pressure-volume diagram. 

In addition, logar ithmetic plots of P, V data provide estimates 
of the combustion delay time, the duration of effective heat 
release, and the ratio of specific heats, y, during the 
isentropic compression and expansion phases. Finally, pressure- 
volume data is required as an input for the thermodynamic models 
used to compute heat release and fuel burning rates. 

Logarithmetic Pressure Volume Diagrams : 

A matrix of comparable test data for engine operation on 
iso-octane and 100-600 is presented in Table 6 and logarithmetic 
P-V diagrams are shown in Figures 40 through 57. When pressure- 
volume data is plotted on a logarithmetic diagram, the isen- 
tropic portions of the compression and expansion process appear 
as straight lines. The slope of the linear segments is equal to 
(-1) .y, where y is the ratio of specific heats. The beginning 
and end of combustion are marked by a departure from and return 
to the straight isentropic compression and expansion lines; 
since the effect of combustion is equivalent to heat being added 



-31- 



wit h a consequent change in y. These points are indicated in 
Table 6. The apparent burning time for the iso-octane and 
100-600 fuels is nearly the same and decreases with increasing 
RPM. The apparent average burning time is 5.2 ms at 1500 RPM, 
3.5 ms at 2000 RPM and 2.7 ms at 2500 RPM. Decreases in 
engine load only slightly decrease the time required for 
burning . 

The combustion delay time can be determined if the start 
of injection is known. The injection is tabulated in Table 6 
and indicated in each figure. In previous work, this delay has 
been attributed to the required jet transit time of the fuel 
from the injector to the stable flame front established at the 

(3) 

spark plug electrodes . However, as shown in Table 6 the fuel 
used also influences the delay time. This indicates droplet 
evaporation rates may also influence the combustion process. 

The small dip in the log P log V diagrams during the injection 
period also indicates the effects of fuel vaporization and this 
dip is more pronounced when the engine is operated on iso-octane 
due to its higher latent heat of vaporization. 

Heat Release Calculations : 

The rate of fuel burning is a basic parameter in most 
engine models and techniques for calculating fuel burning rates 
from engine pressure data have received considerable attention. 
Our understanding of combustion in homogeneous fuel air mixtures 
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is well developed; however, the increased complexity of 

heterogeneous charge engine combustion precludes a strictly 

thermodynamic solution which does not account for mixing in the 

combustion chamber. This statement applies to both stratified 

charge and diesel engines and most of the previous 

work has been done with diesel combustion systems. Until 

recently, efforts to predict heat release rates in a diesel 

engine were highly empirical. Lyn has developed relationships 

to predict heat release rates in an open chamber diesel based 

on the fuel injection rates . Borman and Kreigher have 

developed a thermodynamic model to predict burning rates from 

( 12 ) 

diesel engine pressure data . In the Borman model, the fuel 
was assumed to be homogeneously mixed at each time step. This 
assumption implies very lean equivalence ratios at the start 
of injection which increase to the overall equivalence ratio. 

This is physically inconsistent since mixing considerations of 
the fuel jet imply initially rich combustion followed by pro- 
gressive mixing down to the final lean overall equivalence ratio. 
Still other investigations have assumed micro mixtures of 
burning droplets in which all combustion takes place at stoichio- 

( 13 ) 

metric considerations . The TCCS performance model developed 
by Jain assumed that burning took place at equivalence ratios 
determined by jet mixing and air entrainment and used a specified 
constant equivalence ratio for the burnt products during its 
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control the shape of the burning rate diagram computed from 
pressure-volume data. 

The thermodynamic model used in these experiments to predict 
the cumulative mass of fuel burnt from pressure-volume data is 
an extension of the two zone homogeneous charge combustion model 
as outlined below. The closed system is defined as all air, 
residual gas and fuel vapor in the cylinder prior to ingition. 

The mass conservation equation can be written as: 



v = v/M = xv, + (l-x)v 
b i 



u 



4.1 



and the first law as: 



e = (E -W-Q)/M = xe, + (l-x)e 
o b u 



4.2 



where 



x = charge mass burnt/total charge mass 



= total energy of the charge at time t^ 

M = the total charge mass, air + fuel + residual gas 



Q = the cumulative heat since t 



o 



v = the combustion chamber volume 



W = the work done since t 



o 



e = the average specific energy 
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v = the average specific internal energies 
e^,e^ = the appropriate average specific volumes 
v^,v = the appropriate average specific internal energies 

Subscripts : 

b refer to the burnt zone 
u refer to the unburnt zone 

Further more at a given equivalence ratio $ ; 



v = 
u 



v (P,T ) 
u u 



4.3 



- (P,T ) 
v = v b 

b b 



4.4 



e = e (P,T ) 4.5 

u u u 



e 



b 




cp-V 



4.6 



where 

P = the cylinder pressure 
T^,T = the appropriate average zone temperatures. 

Assuming the unburnt zone undergoes adiabatic quasistatic 
compression and expansion, is calculated from; 



(— ) = r v (P T) - ( — — ) 
op ; s 1 u^l i; dp )r 



4.7 
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Equations 4.1 and 4.2 can be combined to eliminate x, T is 

b 

then found by iterative technique. Once T, is known, x, 

b 

the mass fraction burnt can then be calculated from either 4.1 
or 4.2. 

This model has been expanded by Martin for use in a 

heterogeneous charge engine. The burnt zone is considered as 

burnt products 4- residual gas uniformly mixed at an externally 

designated equivalence ratio <j>. at each time step. The unburnt 

b 

zone is divided into two components; unburnt air + residual gas, 
and the injected but unburnt fuel vapor. Heterogeneous com- 
bustion models require a relationship between the fuel fraction 
burnt and the charge mass burnt since combustion takes place at 
other than the overall equivalence ratios. The following ratios 
can be defined: 

y = unburnt fuel mass/charge mass 
z = fuel mass burnt/total fuel mass 

The equation for specific volume and specific energy of the un- 
burnt zones are then written as: 



v = 

u 



(>v uf + (l-x-y))v ua )/(l-x) 



4.8 



e = 
u 



(y e uf + d-x-y) e ua )/(i-x) 



4.9 
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Equations for y and z are then expressed in dimensionless ratios. 



_ F<j> (1-R ) 
7 1 + F<}> 



F4> b (l-R) 

1+ RF^ + (1-R)F<J> 6 



<)> b (l+F*) 

Z = 

<p (l+RF* + (1-R)F4>J 

b 



4.10 

4.11 



where 

<p, = the burnt zone equivalence ratio at time t 

b 

(p = the average overall equivalence ratio at time t 
F = the stoichiometric fuel-air ratio 
R = the residual fraction 

The computational procedures used to solve for x is the same 

as in the homogeneous charge model; and z can then be found with 

Eauation 4.11; however, particular attention must be paid to the 

change in and <j>_ with time. The overall equivalence ratio only 
b 

varies during the injection process, and its change is directly 
related to the fuel injected in each time step. But it will be 
shown that proper selection of the burnt zone equivalence ratio 
is not straightforward. 

As originally developed by Martin, the modified model as- 
sumed a constant unburnt equivalence ratio equal the final over- 
all equivalence ratio and that fuel burned immediately upon in- 

(4) 

jection . Early efforts to use the model in this form with <j>, 

b 

held constant predicted a slow start of combustion prior to 
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the actual rapid combustion phase, and a final burnt fuel 
fraction greater than one. The model has been refined by in- 
cluding the variations in unburnt equivalence ratio during 
injection, and by assuming that the injected but unburnt fuel can be 
treated as fuel vapor. Procedures for allowing the burnt zone 
equivalence ratio to change with time have also been included 
in the computer program. 

Sensitivity of the thermodynamic model to <j>, is shown in 

b 

Figure 58. These computations assume no mixing, and the burnt 

gas equivalence ratio is held constant. The computations are 

b 

based on pressure-volume data and on overall equivalence 

ratio of i = 0.45 at 1500 RPM. Four values of are shown, the 

b 

overall average equivalence ratio, stoichiometric and two rich 
mixtures. When burning is assumed to take place at the overall 
equivalence ratio $ = 0.45, the fuel fraction burnt does not 
reach 1.0 and the burning time is much larger than predicted 
by the log P - log V diagram. This homogeneous premixed case is 
clearly one limiting example. With either rich or stoichiometric 
burnt gas equivalence ratio, the rapid combustion phase does not 
significantly vary, however, the curves diverge widely and exceed 
a value of unity during the mixing controlled combustion. 

During this period, the mixing rates and burning rates are 
comparable, by definition, and a mixing model is clearly needed 
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to account for the change in <|> due to air entrainment by the 
burning plume if the mass burned is to be correctly related to 
the physical processes in the engine* 

The fuel-air mixtures in the combustion chamber is hetero- 
geneous, and combustion is not limited, a priori to any single 
equivalence ratio at a given time, the burnt and unburned gases 
may continually mix throughout the combustion process. As 

already noted the choice of <{> for each time step during the 

b 

mixing controlled combustion phase is important and an accurate 
entrainment model is required. The air entrainment model pro- 
posed by Blizard and Keck and developed for the TCCS engine by 

(3) 

Jain is used to calculate the mass burned for the same PV 
data used in the sensitivity study as shown in Figure 59. En- 
trainment rates predicted by this model were high and the over- 
all equivalence ratio was reached in 15 CA° , introducing an 
unrealistic dip in the mass fraction burnt curve. As indicated 
by the two previous examples, the calculation of burning rates 
in a stratified charhe engine requires additional development 
to include mixing before plausible results will be obtained. 

We postulate that a model developed to predict NO may serve as 
a tool to aid in untangling the mixing phenomena. The mechanisms 

of NO formation in homogeneous charge engines are relatively 
(14) 

well understood ; and the extension of a homogeneous model to 
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stratified charge engines appears to be plausible. 

PV^ Results 

V 

An examination of the value of PV just before, during and 
after combustion gives a very good idea of the net heat input 
or output to the working fluid due to heat release as a result 

(3) 

of chemical reaction and/or heat loss . Figure 60 is a plot 

of PV Y for the same pressure time data as analyzed in Figures 

58 and 59. The cylinder pressure and volume at the start of 

injection is used as the reference (P V ) . The values of the 

o o 

specific heat ratio before and after combustion were determined 
from Figure 41. The solid line represents heat release at 
constant y^ while the dashed line represents heat release at y^ . 
When the end boundary conditions are applied, namely that the 
process must start as unburnt and end with its maximum coinci- 
dent with the burnt maximum, the lines define a very narrow 
region within which the true heat release curve must fall. The 
dotted line represents the results of the thermodynamic program 
discussed in the last section with the burnt products equivalence 
ratio held at stoichiometric. The rise of this line above the 
peak value evident in the burnt curve is explained by the rise 
in the value of fuel mass fraction burnt in Figure 58 to greater 
than 1.0. 

Several conclusions can be drawn from plots of this nature. 
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The dip in the unburnt curve during the injection period pro- 
vides further evidence of fuel vaporization, a fact also dis- 
cussed in the section covering logarithmetic data plots. If the 
above curves are normalized using the maximum difference in 
burnt and unburnt curves, a line starting as unburnt and changing 
to the burnt curve near TDC closely approximates the actual 

y 

cumulative fuel fraction burnt curve. PV plots with y de- 
termined from log P, log V plots can be used to qualify results 
from a more complete thermodynamic analysis. Note that when 
normalized PV curves are compared to the results of the 
Figure 58, only the fuel fraction burnt curves in which 0.95 < 1.1 

during the rapid combustion phase fall within the defined 
boundary region. 

The method discussed above can, with careful normalization, 
provide a very good estimate of cumulative heat release. These 
estimates can, by comparison with data from detailed thermo- 
dynamic, be used to qualify assumptions made in the necessary 
mixing models. 
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V CONCLUSIONS AND RECOMMENDATIONS 



1. The multifuel capability of the TCCS concept has been 
demonstrated by tests with iso-octane, a wide boiling 
point fuel and methanol. The thermal efficiency of the 
engine is independent of the fuel used. However, proper 
matching of the fuel injection system is required to 
achieve satisfactory exhaust emissions levels. 

2. The limits on engine operating range are determined 
by hydrocarbon and CO emissions. The equivalence ratio 
upper limit is determined by a "smoke limit" near stoichio- 
metric and the lower limit near $ = 0.3 is determined by 
cyclic variations and high hydrocarbon emissions. 

3. Detailed emissions data has been acquired for the 
three test fuels. Emissions with iso-octane and the wide 
boiling point fuel exhibit similar trends and compares 
favorably with previous available data. 

4. Further research is required to explain the effects of 
RPM and different fuel types on indicated mean effective 
pressure at high load conditions with the TCCS system. 
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5. Techniques for obtaining online digital data with a 
small laboratory computer have been demonstrated. The 
pressure-volume data obtained has been shown to be suf- 
ficiently accurate for use in performance models. 

6. It has been shown that accurate pressure-volume data 

can be used to provide a good estimate of cumulative heat 

v 

release through the use of logarithmetic and PV plots. 

7. The problems involved in predicting heat release rates 
have been discussed and the importance of mixing in diesel 
and stratified charge combustion clearly demonstrated. A 
detailed thermodynamic analysis of burning rates will 
require better modelling,- both for the mixing of the fuel jet 
with air before it is entrained in the flame front and for 
the entrainment of air by the burning plume. 

8. It is recommended that a NO prediction model be 

x 

developed for this engine. This model will aid in under- 
standing the role of mixing in stratified charge combustion 
and can be used to qualify mixing models developed for heat 



release calculations. 
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9. Parametric studies involving off design operation are 

required to explain the sensitivity of hydrocarbon emissions 
to small variations in injection and ignition phasing. 
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APPENDIX I 



PERFORMANCE AND EMISSIONS DATA REDUCTION PROGRAM 

An interactive data reduction program was written for the 
Sloan Laboratory PDP 11/10 data analysis facility. The program 



TCCS engine; however, modifications for other single cylinder 
engines are possible. All inputs are requested in the same units 
that are used on the experimental data sheet. Equations used to 
compute air flow rate in grams/second, coolant flow rate, vol- 
umetric efficiency, brake horsepower, and engine emissions require 
clarification. 

Air Flow : 

The equation for the flow rate of air through the ASME square 
edged orifice meter with flange taps is given by the following 
equation (15). 



consists primarily of I/O and is designed specifically for the 




Gy AP 



w = mass flow rate grams /second 



= orifice diameter inches 



K = flow coefficient 



Y = expansion factor 



P^ = static pressure before orifice in in. Hg 
T^ = temperature before orifice °R 
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G = specific gravity of gas 
y = super compressibility factor 
p = pressure drop across orifice in in. H^O 
The computer form of this equation will work with two orifice 
diameters = 0.515 and = 0.71. The program assumes the 
following variable values: 

Y = 1 - 7.3 x 10“ 4 Ap 

r = r r 1+w i 

wet dry L l+1.608w J 

Y « 1 

The maximum and minimum Reynolds number can be written as 

functions of Y and D n with N =0.85, RPM Max = 3000 RPM 

2 v 

Min = 1000, and inlet air temperature = 90°F 

17174.72 „ 68701.3 

Re # = — Re = — 

mm Max 

and appropriate values of K as a function of determined 

0.0366 

K = 0.6152 D 2 
Coolant Flow Rates : 

The equation for coolant flow rates was determined by linear 
regression analysis of calibration data points. 46°F water was 
used and a density correction was applied to obtain the least 
squared curve fit given below, good at 175°F. 
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m = 0.0479 (h) - 0.0081 
r = 0.994 

m = lbm/sec 

h = rotameter height 

r = goodness of fit 



Volumetric Efficiency : 

The engine volumetric efficiency is strictly a function of 
engine dimensions and operating conditions as shown below: 



V 

v 



m( i > 

p i v 

( TjR ) 



m = air flow rate in lbm/min 
N = revolutions/min 
P = inlet air pressure 
V = cylinder volume 
R = specific gas constant 
= inlet air temperature 

after the inclusion of engine geometry and unit conversions 
3 . TDTm (T +460) 

7 ] _ a 1 

v N[0. 0193 P _ - 0.361 P T ] 

atm I 
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Engine Power Output : 

The mean effective pressure and horsepower are determined 
from the following equations (16). 

hp _ me P LAN 
66000 

NAh 



After including specific engine geometry and an overall 
dynamometer constant of K = 6000 



mep = 2.88847Ah 



hp 



NAh 

6000 



where 

c . . . 2 
A = area of engine piston, in. 

L = stroke, ft. 

N = engine RPM 

mep = mean effective pressure, psi.. 

Ah = dynamometer scale height, in Hg. 



Specific Emissions : 

The average exhaust composition is a function of the 
equivalence ratio. A model developed by Stivender ( 8 ) is used 
to determine the exhaust based equivalence ratio as a check of 
the equivalence ratio measured from inlet flow rates; and to 
compute the engine emissions in grams of pollutant/indicated 
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horsepower hour. Required inputs are the indicated specific 
fuel consumption (ISFC) , emissions concentrations on a volume 
basis and the fuel carbon: hydrogen ratio. The model as pre- 
sented does not apply to alcohol fuels and a method presented by 
Spindt 0-7 ) was used for the methanol experiments. 

The model fits one undetermined equilibrium constant for the 
water-gas reduction to direct measurements: 

[HO] [CO] 

K = [co 2 ] [h 2 ] = 3,8 

The combustion reaction for a typical hydrocarbon fuel with air 
may be expressed in the following form: 

CH y + (n) 0 2 + (3.76n)N 2 + (a) CC> 2 + (l-a-6c) CO + 

+(b) H 2 0 + (c) C 6 H 14 + (d) NO - (y/2-b-Tc)H 2 + 

(n-a-(l-a-6c)/2 - b/2 -d/2) 0 o + (3.76 -d/2)N 

Z n Z 

The molecular weight of fuel as it appears in the above 
equation can be written as 

M = 12.01 + 1.008 y 

The molecular weight of air is assumed as 

M = 28.96 
a 
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The air fuel ratio can then be written by a carbon and 
oxygen balance as 



A M 

f = 4 - 76 ir' 

a 



[CO] + [HO] + [NO] 
[C0 2 ] 4- [0 2 1 j- 

[HC] + [CO] + [C0 2 1 



The HC concentration is measured wet. All other pollutant 
concentrations are dry and must be corrected by the following 
relationship 



[ ] 



wet 



- 1 



’dry (1 -' H 2 0,) 



An empirical correlation used to determine the exhaust water 
concentration where the concentrations were wet and K = 3.8, as 
shown below 



c. /[CO,] - [CO] 

'V’ ■ t[c§] 



( 



3 . 8 [ C0 2 ] 



+ 1 ) 



Specific pollutant emissions can be written as 

M [X] 

IS "X" (gr/ihp-hr) - t tHcl + "; C0 ) + CC o 2 ] > ISFC 

where M X" is the species of interest. When the above equation is 

used to indicate C.H . emissions the [”X ,f ] term is [ HC ] / 6 since 

o 14 

[HC] is determined by a count of single carbons. 



- 52 - 



For output consistency fuel consumption is based on the 
observed air flow and the equivalence ratio calculated from the 
exhaust products. The computer program listing is included in 
Appendix III. 
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APPENDIX II 

ONLINE DIGITAL PRESSURE DATA ACQUISITION 

Accurate pressure volume data was required and an online 
digital acquisition system was developed for these experiments. 

The advantages of direct data acquisition include improved accuracy 
and speed. With these routines, a large number of data records can 
be collected and statistically analyzed. Consequently, failure 
of the experimental techniques can be detected by immediate review 
and preliminary analysis of the digital data. The online data 
acquisition system and program described in this appendix is based 
on a Digital Equipment Company (DEC) 11/10 computer with 16 K of 
core memory, a DEC RK -05 random access disk, two teletype 
terminals and the DEC Laboratory Peripheral System (LPS) . The 
DEC LPS consists of an 8 channel multiplexed analog to digital 
converter with variable gain preamplifiers, and internal timing 
clock, and two Schmidt triggers. This system has the capability 
to sample each channel along with its multiplexed pair simul- 
taneously and then perform sequential converstion on the two 
signals. The sample window length is 5 nano-seconds and each 
signal conversion takes 25 yS. Maximum sampling rate on a single 
channel is 45 Hz, with 12 bit conversion. 

A schematic of the data acquisition system is shown in 
Figure 17 . The reference marker pulse at the start of the compression 
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stroke is input to channel 0 with the cylinder pressure input to 
channel 10. Marker pulses every 5 CA° are used as an input to the 
Schmidt trigger. It appears that further system refinements will 
permit 1 CA° sampling increments for a single pair of inputs. 

Accurate pressure records must be matched with accurate 
crankangle records to be of further use and variations in crank- 
shaft angular velocity precludes accurate determination of the 
cylinder volume when the pressure signal is only logged against 
time. Consequently, combustion pressure and crankangle position 
signals make up a data pair and are sychronized with the aid of 
the Schmidt trigger. Each sample interval consists of 144 data 
pairs comprising two engine revolutions. The actual crankangle 
position is not known at the start of a sampling interval and 
the sampled data is reordered at the end of each data sample in- 
terval using a reference signal at 185° before top dead center. 

Two techniques can be proposed for obtaining average engine 
performance. The first involves a complete heat release 
analysis of individual combustion records followed statistical 
averaging of these results. The extended interval between 
sampled data sets and the computational expense required by this 
method precludes its use. Consequently, a second technique in- 
volving the computation of mean cylinder pressure records from 
consecutive data sets was used. These mean records are then 
analyzed to obtain engine performance. If a Gaussian distribution 
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is assumed, the probability density function can be expressed as: 

1 -cv* )2 

p(X) = — — — exp [ x ] 

yj 2 tto 2 a 



The following equations are used to calculate the statistical 
properties of data variations from N records of data: 

N 

Mean (d. ) X. = ^ / X. . 

1 i N 1 3 

J “ i 



X_^ represents the average value of the i the element of the data 
vector X. 

N 



Variance : 



(a 2 ) 






2 



Rearranging these equations, the standard deviation can be 
computed with a single pass of the data: 




This routine permits calculation without an unmanageable number 
of element arrays. However, the standard deviation as calculated 
is not normalized and can only be expressed as a percent of the 
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observed element mean. 

Assuming that we are sampling from a normal population 

it is possible to construct exact confidence intervals for y, 

the true mean, even when a is unknown, by use of the Student - t 

distribution. A 1-a confidence interval for y is expressed as: 

i 



X. 

1 



S. 




V . < x. + t /0 

1 i a/2 




) 



For a large sample size the distribution of can be closely 
approximated as normal and a 1- a confidence interval is: 



S. 

l 




< a . < 

i 



S. 

i 




where 



Z. 

l 



s. - o. 

i l 



(ajyj 2N~ ) 



In addition to the statistical information outlined above, 
the outline data acquisition program integrates the mean pressure 
volume diagram to calculate the indicated mean effective pressure 
and the pumping mean effective pressure by using Simpson’s Rule 
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for non evenly spaced ordinates as shown below: 

N 

W- ^2 [(V j+2 (9) " V j (9) (5P j (9) + 8P j+l (9) " P j+2 (9)] 

j * 1,3,5.* 

Note that cylinder volume is a geometric function of crankangle 
position. 

Log P and Log V vectors are also displayed for cycle eval- 
uation and the program listing containing in stream documentation 
is included in Appendix III. The assembly language commands 
for the sampling subprogram are explained in Ref erence(L8 ) . 
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APPENDIX III 

COMPUTER ANALYSIS PROGRAMS 

Listings for all computer programs and subroutines used 
in these experiments are included in this appendix. The programs 
and subroutines can be grouped in four areas as indicated by 
Tables A-l through A-4 . All programs contain in-stream 
documentation of major equations and computational schemes 
and each subroutine contains a brief section describing its 
purpose, calling sequence, and the definition and dimensions of 
arguments. All programs with the exception of subroutine 
SAMPLE are written in ANS Fortran IV. Subroutine SAMPLE is 



written in DEC Assembly Language, 
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Table A-l 

Summary of Interactive Data Acquisition and Reduction Programs 



Program 


Purpose 


REZLTS 


Calculates basic performance 

and emission^ from experimental data 



ONLINE 


Performs online pressure data ac- 
quisition and calculates mean pressure 
crankangle statistics 



Subroutine 

SAMPLE 


Provides assembly language commands 
used to control analog to digital 
conversion 
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Table A- 2 



Summary of Pressure Crankangle Data 
File Preparation and Control Programs 



Program Purpose 

ANALIZ Prepares pressure crankangle data 

files and basic performance in- 
formation for additional thermo- 
dynamic analysis* 



Subroutine 

RECALL Returns pressure-crankangle data 

from ONLINE for further calculations. 



XPRNT1, 

XPRNT2 



Provides printout of internal 
variables in XCLC2 for each time 
increment . 



STORE 



Provides file storage options 
for output from ANALIZ. 
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Table A- 3 

Summary of TCCS Combustion Analysis Program 



Subroutine 


Purpose 


XCLC2 


Calculates fuel fraction burnt, and 
PVY versus crankangle from pressure- 
crankangle data 



GASVEL 


Calculates the average gas velocities 
at the periphery of the piston cup* 



HEAT 2 


Calculates heat transfer rate in 
TCCS engine during combustion 



PLUME 


Calculates air entrainment rates 
for burning gas plume** 



* Appendix B, Reference (4). 

**Appendix II, Reference (3). 
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Summary of 



Subroutine 

AFTEMP 

CLDPRD 

DERIVS 

HPROD 

TEMP 

TSUBU2 



UPROPZ 



Table A-4 

General Thermodynamic Property Subroutines* 



Purpose 



Calculates adiabatic flame 
temperature from given initial 
state 



Calculates burnt gas properties 
at low (<1000°K) temperatures 



Calculates derivatives of 
properties; for HPROD 



Calculates burnt gas properties 
at temperatures >1100°K 



Calculates T(h,P) for burnt gas 



Calcualtes T(P) from given initial 
state of unburnt gas following an 
isentropic process (assumes no 
fuel vapor) 



Calculates properties of unburnt gas 
(assumes no fuel vapor) 



* Appendixes C and D; Reference (4) 



THIS PROGRAM WILL REDUCE THE EXPERIMENTAL DATA AS IT IS 
REQUESTED. ALL INPUTS WILL BE PROMPTED. 
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TYPE 900 

FORMAT (//////////' ON LINE DATA PROGRAM SUMMARY 



TYPE 901 t IRUNr DAT ( 1 ) r DAT ( 2 > r DAT (3) rRPM 

901 FORMAT ( ' RUN NUMBER ' r 12 r 5X , M r A4 > A4 r 5X , F7 . 2 r 3X , ' RPM 

TYPE 902 

902 FORMAT ( / ' THETA VOLUME/VOLUME MAX MEAN PRESSURE 

1 'STAND DEVIATION LOG(VOL/VOL MAX) LOG(PRESS ATM)') 
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ACCEPT 906 » I TALK 
IF ( I TALK .EG. 1 ) GO TO 500 
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LPSADD-J. 70402 5 VECTOR BUFFER ADDRESS 
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USE SIMPLE ROUTINE FOR LOW TEMPERATURE MIXES 



IF' ( NOTCLD ) GO TO 

CALL CL IT PRD ( P y T y PI 1 1 y DEL r PSI y EN ITTLP y CSUBP t CSUBT * Rl 10 y .DR 1 1 OUT y 
DRUG DP y IER ) 

RETURN 
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GUT PARTIAL DERIVATIVES DY WAY OF A SUBROUTINE CALL HPRl.i.1 360 

I I > I \ ( I I .%/'() 
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Table 1 

Engine Specification s 
DIMENSIONS* 



bore 


3.875 in. 


stroke 


3.875 in. 


connecting rod 


6.625 in. 


clearance volume 


4.570 in.^ 



VALVE TIMINGS 




OPENS 


CLOSES 


inlet valve 10 BTDC (1) 


55 ABDC (2) 


0 (2) 


45 (1) 



exhaust valve 55 BBDC (1) 


10 ATDC (2) 


45 (2) 


0 (1) 


(1) at 0.006 in, valve lift (2) 


valve face flush with head 



defined in Figure ( ) and Appendix B of ( ) 



Table 2 



Injection System 



PUMP 



pump : 
cam: 
plunger : 
reaction valve: 



APE - B Bosch 
6/1 
7mm. 

20 mm^ 



NOZZLE 

nozzle: Roosa-Master XNM 1029 
orifice diameter .023 in. 

orifice L/D 1.0 

nozzle cracking pressure 2000 



needle lift 



.010 in. 
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Table 3 

Summary of Instrumentation 



Air Orifice Inlet 
Air Inlet 
Water Inlet 



Temperatures 
Water Outlet 
Exhaust 
Crankcase Oil 



Bearing Oil 
Fuel Inlet 
Fuel Returns 



Instrument - All Points 

Chromel - Alumel Thermocouple 
Omega DS - 500 Digital Readout 

Resolution 1°F 



Inlet Air 
Crankcase Vacuum 
Oil Pressure 
Exhaust 

Dynamometer Load 
Injection Line 



Combustion Chamber 



Pressures 

Method 

Water Manometer 
Water Manometer 
Panel Gage 
Mercury Manometer 
Mercury Manometer 
Kistler 601 Piezoelectric 
transducer, Kistler 504E 
Charge Amplifier 
Kistler 609A Piezoelectric 
Transducer, Kistler 503D 
Charge Amplifier 



Resolution 
• 1 in . 
.1 in . 

2 PSI 
,1 in. 
.1 in. 

30 PSI* 



.2 PSI** 
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Table 3 (cont) 



Flow Rates 



Air Inlet 
Fuel 

Cooling Water 
Position 
Crankangle 



Method 

ASME Square Edged Orifice 
with water manometers 
Laboratory Scale and Timer 
Rotameter 



Resolution 

0.05 g/sec 
0.01 g/sec 
.2 lbm/sec 



Trump Ross Rotary Pulse Generator 

720 pulses per revolution 

plus Marker .2 CA° 



Fuel Injector 

Needle Lift AVL NH1 - 100 B LDT 1 ym 



Gas Analysis Cart 



Exhaust 
Hydrocarbons 
Nitric Oxides 
Carbon Dioxide 
Carbon Monoxide 
Oxygen 



Scott Model 215 FID HC Analyzer 

TELCO Model 10 A Chemilumenscent NO Analyzer 

Beckman Model 315A NDIR CO^ Analyzer 

Beckman Model 315A NDIR CO Analyzer 

Scott Model 150 Paramagnetic 0^ Analyzer 



* Accuracy is effectively limited by the system transfer 
function and the 565 oscilloscope. 

** See Text 
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Table 4 

Summary of Fuel Properties 



Molecular Wt. 

H:C Ratio 

Specific Gravity 

Boiling Point °F 

Lower Heating Value (Btu/lbm) 

Stoichiometric F/A Ratio 

FIA - 

Aromatics 

Olefins 

Saturates 

Octane No. RON 



Methanol 


Iso-Octane 


Cross-cut 


ch 3 oh 


C 8 H 18 


Distillate 


32 


114 


VL25 


4:0 


2.25 


1.828 


.796 


.692 


o 

00 


149 


211 


(106-648) 


8580 


19080 


18038 


0.155 


0.0665 


0.0692 


% 




29.5 






3.0 






67.5 


106 


100 


76.6 



Cetane No. 



28 ,3 



-137- 



Table 5 





Engine Operating 


Conditions for 


Emission Results 




RPM 




0 

s 


IMEP 


ISFC 


1500 


0.759 


-24 


104.7 


0.384 




0.625 


-23 


93.2 


0.358 




0.446 


-23 


81.9 


0.290 




0.278 


-23 


56.8 


0.284 




0.261 


-23 


53.6 


0.263 




0.172 


-23 


39.1 


0.249 




0.143 


-20 


24.7 


0.331 




0.113 


-23 


17.5 


0.379 


2000 


0.828 


-26 


122.2 


0.353 




0.650 


-26 


108.9 


0.320 




0.475 


-26 


91.3 


0.289 




0.348 


-24 


73.1 


0.269 




0.218 


-23 


44.2 


0.285 




0.161 


-24 


25.9 


0.362 


2500 


0.811 


-29 


130.3 


0.349 




0.639 


-28 


119.3 


0.349 




0.524 


-28 


105.4 


0.286 




0.371 


-28 


82.0 


0.262 




0.249 


-28 


55.5 


0.265 




0.193 


-28 


38.1 


0.306 




0.124 


-28 


21.7 


0,344 



9 = Start of Injection 

s 

Ignition Always Preceeded by 2 CA 
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RPM 




1500 


0.924 




0.762 




0.600 




0.442 




0.426 




0.282 




0.211 




0.129 


2000 


0.927 




0.859 




0.692 




0.472 




0.334 




0.258 




0.180 




0.184 


2500 


1.000 




.850 




.685 




.546 




.367 




.291 




.212 




.218 



Table 5 (con't) 

Cross Cut Distillate 
0 

s 

-22 
-20 
-20 
-19 
-20 
-19 
-19 
-18 
-25 
-23 
-23 
-23 
-23 
-23 
-23 
-20 
-26 
-26 
-25 
-26 
-25 
-26 
-23 
-20 



IMEP 


ISFC 


112.1 


0.432 


105.4 


0.392 


93.0 


0.358 


79.7 


0.299 


76.3 


0.317 


55.5 


0.284 


41.3 


0.290 


22.8 


0.327 


112.9 


0.416 


112.4 


0.390 


98.5 


0.372 


88.4 


0.296 


63.6 


0.303 


51.4 


0.290 


34.9 


0.304 


32.4 


0.338 


122.5 


0.439 


120.7 


0.391 


108.9 


0.353 


95.6 


0.324 


69.3 


0.303 


56.6 


0.300 


39.7 


0.314 


19.4 


0.387 
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Table 5 (cont) 
Methanol 



RPM 




9 

s 


IMEP 


ISFC 


1500 


0.767 


-25 


106.9 


0.878 




0.636 


-28 


101.7 


0.755 




0.486 


-26 


89.0 


0.676 




0.394 


-26 


77.7 


0.637 




0.324 


-26 


66.4 


0.616 




0.261 


-24 


50.8 


0.655 




0.219 


-24 


40.7 


0.697 


2000 


0.811 


-28 


108.9 


0.900 




0.741 


-28 


106.9 


0.838 




0.654 


-28 


104.6 


0.766 




0.591 


-28 


100.5 


0.731 




0.504 


-28 


92.7 


0.682 




0.414 


-28 


84.6 


0.632 




0.364 


-25 


72.5 


0.643 




0.291 


-25 


58.6 


0.640 




0.221 


-25 


26.9 


1.077 


2500 


0.696 


-35 


106.3 


0.824 




0.641 


-35 


102.3 


0.795 
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Table 6 

Matrix of Averaged Pressure Crankangle 
Iso-Octane 



RPM 






IMEP 

m 


IMEP. 

l 


e. 

IS 


0. 

le 


9 bs 


9 be 


1500 


0.77 


0.76 


104.7 


106.6 


-24 


9 


-12 


35 




0.46 


0.45 


81.9 


82.3 


-23 


0 


-13 


35 




0.29 


0.28 


56.8 


60.0 


-23 


-3 


-12 


24 


2000 


0.83 


0.83 


122.2 


118.3 


-25 


13 


-12 


28 




0.49 


0.48 


91.3 


88.5 


-24 


3 


-9 


24 




0.22 


0.22 


44.2 


44.3 


-23 


-5 


-7 


23 


2500 


0.78 


0.81 


130.3 


122.4 


-29 


15 


-8 


30 




0.52 


0.52 


105.4 


99.6 


-28 


7 


-9 


30 




0.25 


0.25 


55.4 


47.4 


-28 


-4 


-9 


29 



Subscript Nomenclature 

0 = Observed from fuel air flow 

e = Measured from exhaust gas composition 
M = Dynamometer measurement 

1 = Integrated average pressure data 
is = Injection start 

ie = Injection end 

bs = Start of heat release from Log P Log V plots 

be = End of effective heat release from Log P Log plots 
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Table 6 (cont) 
Cross Cut Distillate 



RPM 




*e 


IMEP 

M 


IMEP. 

l 


9. 

is 


0 . 

le 


9 bs 


V 


1500 


0.76 


0.78 


105.4 


99.2 


-20 


10 


-11 


35 




0.43 


0.43 


76.3 


73.3 


-20 


4 


-13 


36 




0.27 


0.27 


51.4 


47.2 


-20 


-5 


-11 


35 


2000 


0.87 


0.85 


121.0 


114.6 


-24 


15 


-14 


29 




0.49 


0.47 


88.4 


84.9 


-23 


-2 


-14 


26 




0.26 


0.26 


51.4 


46.5 


-23 


-12 


-15 


31 


2500 


0.83 


0.82 


121.3 


124.8 


-26 


15 


-14 


31 




0.52 


0.50 


95.0 


100.1 


-26 


0 


-14 


24 




0.29 


0.28 


56.1 


59.7 


-26 


-9 


-13 


26 



//////////////// 
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Fig, 1 -TEXACO CONTROLLED COMBUSTION 
SYSTEM (SCHEMATIC) 



CRACKING PRESSURE ADJUSTMENT SCREW 
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ROOSA-MASTER INJECTION NOZZLE 
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SPARK PLUG 
AXIS 




THE PLUME 
INITIAL STAGE 

COMPLETE 

TORUS 



INJECTOR AXIS 



YZZZ/ZZZZZZZZZZZA 
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SCHEMATIC OF THE COMBUSTION PROCESS 

Fig, 3 
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FIGURE 5 PHOTOGRAPH OF TEST INSTALLATION 
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Fig. 6a Single Cylinder Engine Assembly 



- 148 - 




Fig. 6b Single Cylinder Test Assembly 



- 149 - 




Fig. 7 Dynamometer Hydraulic Scale 
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RESERVOIR 




Fig. 8 cooling system 
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PUMP 



Fig. 9 LUBRICATION SYSTEM 



ATMOSPHERE 
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FLOW PRESSURE NEEDLE LIFT 

TRANSDUCER TRANSDUCER INDICATOR 
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Fig. 12 INJECTOR LIFT TRANSDUCER ASSEMBLY 
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2 . FUEL INJECTOR NEEDLE LIFT 

3 . FUEL INJECTION LINE PRESSURE 

4 . CYLINDER PRESSURE 

5 . ■ crankangle INDICATOR, 1/5 CA° 

6 . 18 5 CA° BTDC REFERENCE SIGNAL 

7 . balance pressure indicator 



Fig.13 -TYPICAL OSCILLOSCOPE TRACE 



LOG PRESSURE (ATM) 
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Fig. 14 1500 RPft Motoring Log.P vs Log V 



LOG PRESSURE ( ATM ) 
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FiG. 15 2000 RPfl Motoring Log P vs Log V 



- 158 - 




LOG(VOLUME /VOLUME MAX) 



Fig, 16 2500 RPM Motoring Log P vs Log V 



- 159 - 




Fig, 17 Schematic of Data Acquisition System 



INDICATED MEAN EFFECTIVE PRESSURE -PSI 
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Fig. 18 Indicated Mean Effective Pressure vs 
Equivalence Ratio, for Iso-octane 



INDICATED MEAN EFFECTIVE PRESSURE -PSI 
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1 




Fig. 19 Indicated Mean Effective Pressure vs 
Equivalence Ratio, for 100-600 Fuel 



- 162 - 




Fig. 20 Indicated Specific Fuel Consumption vs 
Equivalence Ratio, for Iso-octane 



INDICATED SPECIFIC FUEL CONSUMPTION LB/IHP-HR 
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Fig. 21 Indicated Specific Fuel Consumption vs 
Equivalence Ratio, for 100-600 Fuel 
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Fig. 22 Indicated Thermal and Volemetric Efficiency 
vs Equivalence Ratio for Iso-octane 
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Fig. 23 Indicated Thermal and Volumetric Efficiency 
vs Equivalence Ratio for 100-600 Fuel 
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Fig. 24 Exhaust Temperature vs Equivalence Ratio 
for Iso-octane 
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Fig. 25 Exhaust Temperature vs Equivalence ratio for 
100-600 Fuel 



CO EMISSIONS GR/IHP-HR 
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Fig. 26 Carbon Monoxide Emissions vs Equivalence 
Ratio for Iso-octane 
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100-600 D 




Fig, 27 Carbon Monoxide Emissions vs Equivalence 
Ratio for 100-600 Fuel 
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ISO-OCTANE 
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Fig. 28 Hydrocarbon Emissions vs Equivalence 
Ratio for Iso-octane 



HC EMISSIONS (C 6 H, 4 ) GR/IHP-HR 
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FiG. 29 Hydrocarbon Emissions vs Equivalence 
Ratio for 100-600 Fuel 
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Fig. 30 Nitric Oxide Emissions vs Equivalence 
Ratio for Iso-octane 



N0 X EMISSIONS (N0 2 ) - GR/I HP-HR 
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FiG. 31 Nitric Oxide Emissions vs Equivalence 
Ratio for 100-600 Fuel 
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Fig. 32 Indicated Mean Effective Pressure vs 
Equivalence Ratio for Methanol 
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Fig. 33 Indicated Specific Fuel Consumption for 
Methanol 
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Fig. 34 Indicated Thermal and Volumetric 

Efficiency vs Equivalence Ratio for 
Methanol 
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Fig. 35 Exhaust Temperature vs Equivalence 
Ratio for Methanol 
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Fig. 36 Carbon Monoxide Emissions vs Equivalence 
Ratio for Methanol 
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Fig, 37 Hydrocarbon Emissions vs Equivalence 
Ratio for Methanol 
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Fig. 38 n itric Oxide Emissions vs Equivalence 
Ratio for Methanol 
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Fig. Volumetric Efficiency vs Engine RPM 




Fig. 39 Friction Mean Effective Pressure 
vs Engine RPM 
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Fig. 40 Log P vs Log V for Iso-octane, 
* = 0.7G, 2500 RPfl 



LOG PRESSURE (ATM) 
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LOG (VOLUME / VOLUME MAX) 



Fig. 41 Log P vs Log V for Iso-octane, 
+ = 0.45, 1500 RPM 



LOG PRESSURE (ATM) 
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Fig, 42 Log P vs Log V for Iso-octane, 
* = 0,28, 1500 RPM 



LOG PRESSURE (ATM) 
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Fig, 43 Log P vs Log V for Iso-octane, 
j = 0.83, 2000 RPM 
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Fig. 44 Log P. vs Log V for Iso-octane, 
♦ = 0.48, 2000 RPfl 
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- 187 - 




.2 



2000 RPM ISO-OCTANE $ 0.22 




-1.0 -.9 -.8 -.7 

LOG( VOLUME / VOLUME MAX) 



Fig. 45 Log P vs Log V for Iso-octane, 
* = 0.22, 2000 RPM 



LOG PRESSURE (ATM) 
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Fig. 46 Log P vs Log V for Iso-octane, 
♦ = 0.81, 2500 RPN 



LOG PRESSURE (ATM) 
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Fig. 47 Log P vs Log V for Iso-octane, 
♦ = 0.52, RPM 



LOG PRESSURE (ATM) 
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FiG. 48 Log P vs Log V for Iso-octane, 
co = 0.25, 2500 RPM 



LOG PRESSURE (ATM) 
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Fig. 49 Log P vs Log V for 100-600, 
♦ = 0.77, 1500 RPM 



LOG PRESSURE (ATM) 
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Fig, 50 Log P vs Log V for 100-600, 
<0 = 0.43, 1500 RPM 



LOG PRESSURE (ATM) 
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FiG. 51 Log P vs Log V for 100-600, 
♦ = 0.27, 1500 R ?n 



LOG PRESSURE (ATM) 
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Fig, 52 Log P vs Log V for 100-600, 
♦ = 0.85, 2000 RPM 



LOG PRESSURE (ATM) 



- 195 - 




-1.1 -1.0 -.9 -.8 -.7 -.6 -5 

LOG( VOLUME / VOLUME MAX) 



Fig. 53 Log P vs Log V for 100-600, 
♦ = 0.47, RPM 
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Fig, 54 Log P vs Log V for 100-600, 
♦ = 0.26, 2000 RPM 



LOG PRESSURE (ATM) 
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FiG, 55 Log P vs Log V for 100-600, 
♦ = 0.82, 2500 RPfl 



LOG PRESSURE (ATM) 
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LOG (VOLUME / VOLUME MAX) 

Fig. 56 Log P vs Log V for 100-600, 
* = 0.50, 2500 P.Pfl 



LOG PRESSURE (ATM) 
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FiG. 57 Log P vs Log V for 100-600, 
♦ = 0.28, 2500 RPM 
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Fig. 58 Fuel Mass Fraction Burnt vs Crankangle for Fixed 4* 
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Fig. 59 Fuel Mass Fraction Burnt vs Crankangle for Variable 



= f C m b) 



-202- 




<d 

a 



oo 

uu 

LjU 

CD 

LjU 

Q 



<d 



crz 

c_> 



co 

LU 



< 

> 

< 

z: 

z: 

< 

CD 



CK 

Z) 

OQ 



Q 

z: 

< 

H 

z: 

cr 

ZD 

PQ 



LU 

-J 

CO 

z: 

< 

z: 

< 

cc 

CJ 

CO 

> 

>- 

> 

a_ 

o 

CO 

CD 



tAarsh j e ta>' ed P erf °of 
& aeta ‘-on ° T 

nre comp 3 "' » s *, n 

the; “ eng' 06 * 

charge e g 



-****?% 






M)0’ 



Thesis 
M3 547 



1 3 3 1 

M a rsh 

A detai led perfor- 
mance comparison o £ 
disti 1 late fuels in 
they exaco stratified 
charge engine. 



thesM3547 

A detailed performance comparsion of dis 




3 2768 002 12794 6 



DUDLEY KNOX LIBRARY 



